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Abstract 

The semiclassical Euchdean path integral method is apphed to compute the low temperature 
quantum decay rate for a particle placed in the metastable minimum of a cubic potential in a finite 
time theory. The classical path, which makes a saddle for the action, is derived in terms of Jacobian 
elliptic functions whose periodicity establishes the one-to-one correspondence between energy of the 
classical motion and temperature (inverse imaginary time) of the system. The quantum fluctuation 
contribution has been computed through the theory of the functional determinants for periodic 
boundary conditions. The decay rate shows a peculiar temperature dependence mainly due to the 
softening of the low lying quantum fluctuation eigenvalues. The latter are determined by solving 
the Lame equation which governs the fluctuation spectrum around the time dependent classical 
bounce. 
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1. Introduction 



Quantum tunneling from a metastable state through a potential barrier is a fundamental 
nonlinear phenomenon occuring in many branches of the physical sciences 
^, qI- The standard semiclassical approach to metastability has been formulated since long 
by Langer and Coleman 11 1 in the fields of statistical and nuclear physics respectively. 
The main idea underlying such approach consists in selecting a classical background which 
solves the Euler-Lagrange equation and makes a saddle point for the action. Around the 
background, the quantum fluctuations are treated in quadratic approximation and their 
spectrum is obtained by solving a Schrodinger like stability equation whose potential is 
given by the second spatial derivative of the metastable potential. The semiclassical method 
finds a concise and powerful description in the Euclidean path integral formalism 
in which the time for the bounce to perform a full excursion (inside the classically allowed 
region) is a measure of the inverse temperature of the system. In the standard treatments 
of metastability 10, lUl, it is assumed that such time is infinite and therefore the decay 
rate formula holds, strictly speaking, only at T = 0. For applications to specific systems 
however, a precise knowledge of the decay rate at finite T (but within the quantum tunneling 
regime) may be of great interest. To this purpose one has to build the finite time theory 
of metastability for specific nonlinear potentials, setting the crossover temperature between 
(low T) quantum and (high T) activated regimes and find the shape of the decay rate 
when the crossover is approached from below. Focussing on a widely investigated model in 
nonlinear science, a particle in the one dimensional cubic potential, I present in Section 2 
the finite time solution of the Euler-Lagrange equation in terms of the powerful Jacobian 
elliptic functions formalism [14]: this generalizes the well known infinite time bounce which is 
recovered asymptotically. I emphasize that the system is taken as non dissipative therefore 
the temperature should not be viewed here as a property of the heat bath 15, 
but rather as a measure of the system size along the time axis. Section 3 is devoted to the 
computation of the classical action. Section 4 describes the method of the semiclassical path 
integral and presents the calculation of the overall quantum fluctuation contribution through 
the theory of the functional determinants. Section 5 solves the stability equation for the 
periodic potential defined by the classical background. This permits to obtain analytically 
the lowest quantum fluctuation eigenvalues as a function of the finite time/temperature. It 
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is shown in Section 6 that the softening of such eigenvalues close to the crossover largely 
determines the peculiar shape of the decay rate and its deviation from the prediction of the 
standard zero T theory. The conclusions are drawn in Section 7. 

2. Cubic Potential Model 

To begin, consider a particle of mass M in the one dimensional cubic anharmonic poten- 
tial: 

V{x) = - Ix' , (1) 

plotted in Fig. [TJ^a) for fiw = 20meV and M = lO^rrie, me being the electron mass. Say 
a the position of the top of barrier whose height is V{a) = 70^/6 with 7 = Mu'^/a. Let's 
take throughout the paper, a = lA. At x = the particle is in a local minimum from 
which it cannot escape classically. Thus, in the real time formalism, the classical equation 
of motion admits only the trivial solution Xd = 0. Physically however such local minimum 
is metastable as quantum fluctuations allow the particle to explore the abyss at a; > 3a/2. 
In fact, a non trivial classical solution can be found in the Euclidean space. Performing a 
Wick rotation from the real to the imaginary time, t —iT, the equation of motion reads: 



Mxe/(r) = V'{x,i{t)) , (2) 

where V means derivative with respect to Xd- In the spirit of the semiclassical method, 
the particle path x(r) has been split in the sum of a classical and a quantum component, 
x(r) = Xd{T) + ri{r). The Wick rotation is equivalent to turn the potential upside down 
with respect to the real time as shown in Fig. [Dj^b). Now it is clear that the classical motion 
can take place in the reversed potential: precisely, the particle moves back and forth between 
the turning points Xi and X2 in which the particle velocity vanishes. 

Integrating Eq. ([2]), one gets: 

^xl{T)-V{xd{T)) = E , (3) 
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with the constant E representing the classical energy. Defining: 



Xd{T) = 

6 a 
AE 



Eq. ([3]) is easily integrated to yield: 



(4) 



1 [^-^^^ dx 
r - To = ±- / = , (5) 

where tq is the center of motion between the turning points. The boundary conditions 
on the classical motion define a physical picture in which the particle starts from X2 at the 
time r = —L/2, reaches xi at r = tq and returns to the initial position at r = L/2. Then, 
Eq. ([5]) has a time reversal invariant solution whose period L is finite and dependent on E. 

Looking at Fig. [H^b), one sees that the amplitude xi —X2 attains the largest value for the 
E = Q motion while X2 and x^ coincide. For E < Q motions, x^ is negative. The turning 
points Xi, X2, x^ are given by the zeros of the equation — + + k = (x = 2a;/(3a)) 
which admits three real solutions for n G [—4/27,0], that is for E G [— \^(a),0]. After some 
algebra I find: 



Xi= 2 + 2 cos(i9) 

X2= i + ^cos(^-27r/3) 

X3 = ^ + ^ cos(^ - 47r/3) 

1 ,21k . 

v = -arccos(— - — h Ij . (6) 

At the bounds of the energy range, Eq. ([6]) yields: 



E = ^ Xi = 1; X2 = X3 = 

E = -V(a) ^ xi = X2 = 2/3; xs = -1/3 



(7) 
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Thus, at the sphaleron energy Egph = \E\ = V{a), the amphtude of the finite time 



solution has to shrink into a point. Let's find the genera 
center of motion at XciiTo) = Xi ci-nd using the result [l8l |: 



» dx 2F(A,p) 



solution of Eq. ([5]) by pinning the 



A = arcsin I 



Xi - Xcijr ) 
Xi - X2 



P-J^^. (8) 

V xi - X3 

where F{X,p) is the elliptic integral of the first kind with amplitude A and modulus p. 
Then, through Eqs. (j4]), ([5]), ([8]), 1 derive the bounce solution of the finite time theory: 



3ci 

Xdir) = Yixicn'^i^^P) + X2sn'^iTu,p)] 



^= VXl - X3 - To) , 



(9) 



sn{w,p) and cn{zu,p) are the sine- and cosine- amplitudes respectively [l^. The modulus 
p keeps tracks of the classical mechanics through the second of Eq. and Eq. ([6]). 
At E = 0, p= 1, the bounce of the infinite time theory is recovered: 



Xci{r) = Ycn^(ti7, 1)= ysec/i^(^(r - To)) 



(10) 



At Egph, p = 0, the bounce solution is (as expected) a point-like object set at the bottom 
of the valley in the reversed potential: Xci{t) = a. Thus, Eq. ((91) defines the transition state 
which is a saddle for the action below the sphaleron. Computation of Eq. ([9]), shows that the 
bounce amplitude contracts by increasing the energy over potential height ratio (in absolute 
value). As the bounce is a combination of squared Jacobi elliptic functions, its period is 
2K{p) with K{p) = F{t[/2,p) being the complete elliptic integral of the first kind fl^ . 
Hence, from Eq. ([9]), I get: 
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v/^r^^L= i^b) , (11) 

which estabhshes the relation between the oscillation period and the classical energy 
embedded in the turning points. As stated in the Introduction, one can map the imaginary 
time onto the temperature axis, L = h/lKsT*), where T* is the temperature at which 
the particle makes the excursion to and from the edge of the abyss for a given E. Then, 
only periodic bounces whose period is proportional to the inverse temperature determine 
the decay rate and the finite time theory can be viewed as a finite T* theory. From Eq. (fTTj) 
I get: 



4 K{p) ^ ' 

Eq. (fT2!) is plotted in Fig. [2] on a linear scale. Approaching E = 0, T* consistently drops 
to zero while the value at Esph defines the transition temperature T* between quantum and 
activated regimes. 

Analytically, at the sphaleron, Eq. f[T^ yields 



KbT* = — 



2^ ' 



(13) 



which represents the upper bound for the occurence of quantum tunneling and precisely 
sets the Goldanskii criterion 19|, |20] for a cubic anharmonic potential. Taking hu = 20meV, 
I get T* = 36.938i^. The following calculations are carried out in the low temperature range 
up to T*. 



3. Classical Action 



The classical action A[xc/] for the bounce in the finite time theory can be computed in 
terms of the path velocity Xd{T) by the relations: 
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A[xci] = MN-^ - E ■ L{E) 
Jo 

^ = -^{Xi - X2)VXi - X3 , (14) 
where dn{-cu,p) is the delta- amphtude jl^ . 

Computation of Eq. f|T^ requires knowledge of L{E) through Eqs. (jlj), and (fTTj) . 
In the ii^ — > hmit, from Eq. (IT^ . I get the result 



' ^^^^ 

which serves as testbench for the computational method. The dependence of the classical 
action on 1/7^ reflects the well known fact that metastable systems are non perturbative 
and provides the fundamental motivation for the semiclassical treatment. Eq. f[T^ permits 
to set the potential parameters such as the condition A[xc/] > h holds and the semiclassical 
method is thus justified. As M and a have been taken constant, A[xci] oc u in the E —* 
limit. 

The classical action and the squared norm of the path velocity (times M) are displayed in 
Fig. [3]^a) and Fig. [3]^b) respectively. While, at low T*, the two plots are essentially identical 
the role of the term E ■ L{E) in Eq. fll4p becomes more significant at increasing T*. At 
T*, N~'^ vanishes whereas v4[a;c«] is finite. Note that v4[a;ci] decreases smoothly versus T* 
confirming that the transition to the activated regime at T* is of second order as suggested 



long ago 



20 



2l|. 



In general, the criterion to establish the order of the transitions in periodic tunneling 
systems has been formulated by Chudnovsky 22| through the behavior of the oscillation 
period L{E): i) a. monotonic L{E) below the sphaleron implies A[xc«] < Aq for T* < T*, 
with Aq being the thermal action given by Aq = hV{a)/ KbT* . At T* = T*^ both conditions 
^[a^c/] = ^0 and dA[xci\/dT = dA^/dT are fulfilled hence the crossover from the quantum 
to the thermal regime is expected to be smooth; ii) on the other hand, a nonmonotonic 
behavior of L{E) would indicate a sharp transition. As it can deduced from Fig. [2], L{E) 
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increases versus E G [—V{a), 0], thus the case i) apphes to the cubic potential in Eq. ([T]) and, 
consistently, the action is convex upwards versus T*. At the sphaleron, I find numerically: 
L{Esph)/h = 0.314meV^~^. Note, from Eq. fllip . that such value corresponds to 2TT/hw 
(being K{p = 0) = vr/2) and this proves the correctness of the computation. In fact at Egph 
the bounce is a point, that is a static solution of Eq. Q but, near Egph, the periodic path 
is the sum of the sphaleron and a fluctuation with negative eigenvalue e^i whose period 
tends to L{Esph) = 27r/ \/\e-i\ 



2J]. Then one infers that, for \E\ —>■ Esph, the ground 



state eigenvalue has to behave as: e^i —uJ^. This key point will be further investigated 
in Section 5. 



4. Semiclassical Euclidean Path Integral 

This Section presents the calculation of the space-time Euclidean path integral between 
the positions xi and x/ connected in the time L. In the semiclassical model and treating 
the quantum fluctuations in quadratic approximation, the path integral reads: 

< Xj\xi >L= exp 
AfH = J ^^/^(y^'(^) + ^V"{x^,{T))v'{r)j 

V"(x,Kr))= Mc.2(l-^XeKr)) . (16) 

Thus, to get the quantum fluctuation action v4j[?7], one has to solve a second order 
differential problem which, after partial integration in the flrst term, is formulated as follows: 



Ajxci] 

h 



Drjexp 



h 



6 = -dl + V"{x,i{t))/M 

oo 
n=-l 

where the e„ are the quantum fluctuation eigenvalues while the coefficients ^„ of the 
series expansion in ortonormal components ^^^(t) define the measure of the fiuctuation paths 
integration in Eq. ( fT6|) : 
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Dr^=^W / , (18) 

K depends only on the functional integral measure. 

First, observe from Eq. ([2]) that Xci{t) satisfies the homogeneous equation associated 
to the second order Schrodinger-like differential operator Or]n{T) = 0. This is a general 
consequence of the r-translational invariance of the system. Hence, Xci{t) is proportional to 
the ortonormal eigenmode tiqIt), (?7o(t) = Nxcl{t)) with £0=0- The latter however cannot 
be the ground state as the bounce solution (Eq. ([9])) is non monotonic and Xci{t) has one 
node along the time axis within the period L. This implies that the quantum fluctuation 



spectrum has one negative eigenvalue corresponding to the ground state [25|. Here lies the 
origin of metastability. Second, from Eq. ( JT^ . note that for any two points roi, -072 such that 
W2 = zui ± 2K{p), Xci{^2) = icii'^i)- The important consequence is that the fluctuation 
eigenmodes obey periodic boundary conditions (PBC). 

As Xi and Xf defined in Eq. (fT6|) coincide for the periodic bounce, Eq. (fT6|) represents 
the single bounce contribution Zi to the total partition function Zt- In fact, the latter 
also contains the effects of all multiple (non interacting) excursions to and from the abyss 
which is equivalent to sum over an infinite number of single bounce contributions like Zi. 
Moreover, also the static solution of Eq. ([2]), Xc/ = contributes to Z^ by the harmonic 
partition function Z^ which can be easily determined using the same measure in Eq. (ITSl) . 
Summing up, Zx is given by 



Zt = Zhexp{Zi/Zh) 

Z^ = K|Det[/l]|-l/2 ^ (ig) 

Det\}i\ (h = —d'^ + cu^) is the harmonic fluctuation determinant. Being the decay rate F 
proportional to the imaginary exponential argument through the Feynman-Kac formula , 
it follows that there is no need to determine K as it cancels out in the ratio Zi/Zh.. 

Supposed to have solved Eq. ( |T7I) . the quantum fluctuation term in Eq. (|T6l) can be 
worked out by carrying out Gaussian path integrals. Formally one gets: 
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Dt] exp 



AM 

h 



^■Det 



O 



-1/2 



Det[d] = 



(20) 



n=-l 



The evaluation of Eq. fl20l) is carried out through the two following steps. 



A. Zero Mode 



The Gaussian approximation leading to Eq. fl2Ul) is broken by the Goldstone mode arising 
from the fact that tq, the center of the bounce, can be located arbitrarily inside L. The 
technique to overcome the obstacle is well known [20]: the divergent integral over the coor- 
dinate d% associated to the zero mode in the measure Dt] is transformed into a dr^ integral. 
Accordingly the eigenvalue = is extracted from Det[0] and its contribution to Eq. (!20ll 
is replaced as follows 



27lh 

(21) 



To be rigorous, this replacement holds in the approximation of quadratic fluctuations 26] 



while higher order terms may be significant around the crossover. It is also worth noticing 
that Eq. (12T|) is often encountered in the form (eo)^^^^ \/ ^2w'h However the latter is 
correct only in the low T limit where A[xci] equals MN^"^ (as made clear by Fig. [3]) while, 
approaching T*, the difference between the two objects gets large. This fact is crucial in 
establishing the behavior of the decay rate at the crossover as shown in Section 6. 

Thus, handled the zero mode, I turn to the evaluation of the regularized determinant 
Det^p] defined by Det[d] = Eq ■ Det^[d]. 



B. Regularized Fluctuation Determinant 



The calculation of Det [O] is based on the theory of functional determinants for second 
order differential operators which was first developed for Dirichlet boundary conditions 



and then extended to general operators and boundary conditions in several ways 
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29 
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As a fundamental feature, to evaluate Det^\0] one has to know only the classical path 
which makes the action stationary. As shown above the path velocity obeys PBC for any 
two points tui, 'DJ2 separated by the period 2K{p). The latter corresponds to the oscillation 
period L along the r-axis. It can be easily checked that also the path acceleration fulfills 
the PBC. Then, the regularized determinant is given by [29(: 



< /ol/o > (/i(t^2) - /i(roi)) 



/o(ti7i)iy(/o,/i) 

where /o,/i are two independent solutions of the homogeneous equation: 0?7„(r) 
/o is obviously Xd while fi can be taken as: 



(22) 
= 0. 



fi 



dq 



q = p 



(23) 



iy(/o, /i) is their Wronskian and < /o|/o >= A^"^ is given by Eq. f|T4|) . 
The Wronskian, being constant along r, can be calculated in any convenient point. Let's 
take To as /o(to) = 0. Then: 



WifoJi 
9 



-/o(ro)/i(ro) 



^a^uj^iXi - X2)iXi - X3)^ 



Working out the calculation, Det^[0] in Eq. (122!) transforms into: 



(24) 



1-/ , 



E{n/2,p)-p'K{p) 

p2 



< /ol/o > 
WifoJi) 



(25) 



which can be directly computed using Eqs. (HM . ([21]). E{tt/2,p) is the complete elliptic 
integral of the second kind [14 1. 



It is however known in the theory of functional determinants 26|, |27l] that only ratios 
of determinants are meaningful in value and sign, such ratios arising naturally in the path 
integral method as it has been pointed out above. In fact, Det^[0] would diverge in the 
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T* ^ limit due to the fact that the determinant is the product over an infinite number 
of eigenvalues with magnitude greater than one. Consistently with Eq. f|T9|) . Det^[0] has 



to be normalized over Det[h] which, in the case of PBC, is: Det[h] = — 4 sinh^(ct;L/2) 29|. 
The normalization cancels the exponential divergence and makes the ratio finite. 
Then, observing that for E ( T* — 0): 



TO O 
r I r 62 

< /o /o >^ -a ^ 

5 

Kip) ^ ln(4/v/l-p2) , (26) 
from Eq. fl25l) . I finally get the finite ratio 

Det[h] GOco^ ■ ^ ^ 

The dimensionality [i^"^] correctly accounts for the fact that one eigenvalue has been 
extracted from Det^[0]. From the computation of Eq. (125|) . the following informations can 
be extracted: 

1) the T* — > limit given by Eq. ( 1271) is in fact an excellent estimate up to T* ~ T* /2 
whereas a strong deviation is found at larger T* up to ~ T*. 

2) The ratio Det^[0]/Det[h] is negative for any T* and this sign has physical meaning as 
it is due precisely to the negative ground state eigenvalue e_i of the fiuctuation spectrum. 
As Zi/Zh contributes to the partition function by the square root of the fiuctuation (inverse) 
determinants ratio it follows that such contribution is purely imaginary. 

Moreover, close to the sphaleron, T* ~ T*, < /o|/o > oc and Vr(/o, /i) oc p^, thus 
Det^[0] tends to zero as Det^[0] oc p"^. The fact that Det^[0] vanishes at the crossover 
causes the divergence of the inverse ratio. This divergence may look surprising since the 
zero mode had been extracted from the determinant: physically one realizes that there must 
be a quantum fiuctuation mode (ei) which softens by increasing T* and ultimately vanishes 
at the sphaleron. 

To understand in detail the key effects of the low lying fiuctuation eigenvalues ei and 
one has to determine them analytically by solving the stability equation in Eq. ([1] 
This is done in the next Section. 
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5. Lame Equation 



Take Eq. (fT71) with the second derivative of the potential given in Eq. ( IT6|) . Using Eq. ([9]) 
and working out the algebra, I get the stability equation which governs the fluctuation 
spectrum around the classical background: 



_ 4(1 -3xi) 4£„ 



/(/ + 1)^12. (28) 



This is the Lame equation in the Jacobian form for the case / = 3 [31|. For a given I 
and p, Eq. (|28i) yields periodic solutions (which can be expanded in infinite series) for an 
infinite sequence of characteristic An values. The continuum of the fluctuation spectrum 
stems from this sequence. However, being / positive and integer, the first 2/ + 1 solutions 
of Eq. (!28|) are not infinite series but polynomials in the Jacobi elliptic functions with real 
period 2K{p) or 4K{p). Being the period of the potential, 2K{p) plays the role of a lattice 
constant. 

Then, Eq. ( l28l) admits seven polynomial solutions with eigenvalues An,n G from 
which the corresponding are derived. However not all the En are good fluctuation eigen- 
values. In fact four out of seven have to be discarded as their eigenf unctions do not fulfill 
the PBC required for the fluctuation components: ?7„(ci7i) = r7„(ti7i =F 2K{p)). Thus, the 
three good eigenmodes and relative eigenvalues in polynomial form are: 
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rjo oc snl'OJ , p)cn{w , p)dn{tz> , p) 
eo= 



rji oc [sn vj 
El = uo^[ ai 



?7_i oc {sn zu — p 



Ai = -(2 + 5p2) -2v/V 



a^Ai 

2U/2 



sn w + 



P^ + Ai 



sn w + 



p2 + A-i 



p2 + 1 



A-i = -(2 + 5p2) + 2^/Ap^-p^ + 1 

^2 = 



1-3x1 

Xi - X3 



4 



(29) 



The plots of El and e^i versus the energy over potential height ratio are reported on 
Fig. 111(a) and Fig. Hl^b) respectively. Note that: 

i) Eq = is the zero mode eigenvalue correctly recovered through the stability equation. 

ii) El lies in the continuum and, as it can be easily deduced from Eq. fl29l) . it drops to 
zero close to the sphaleron as Ei oc p^: this is precisely the behavior previously envisaged 
for Det^[0]. Hence, Ei is the soft mode driving the enhancement in the decay rate below 
the sphaleron which is discussed in the next Section. Observe that, for p ^ 0, K{p) ~ 
tt/2 + np'^/8. Then, by Eq. ( |T2l) . close to the crossover: Ei oc T* — T*. At T*, E\ and Eq 
merge consistently with the double degeneracy of the corresponding eigenmodes above the 
crossover. 

iii) E-\ is the negative eigenvalue responsible for metastability. also softens (in 
absolute value) with respect to the value e_\ = — 5ct;^/4 found at = 0. Interestingly, 
along the temperature scale, the substantial reduction starts up at T* ~ T*/2, that is in 
the same range at which the classical properties deviate from the predictions of the infinite 
time theory. Finally, at the sphaleron, from Eq. (1291) I get E^i = —uj^ thus confirming the 
prediction made at the end of Section 3. This completes the analysis of the soft eigenvalues 
which ultimately govern the quantum fluctuation spectrum. 
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6. Decay Rate 



The decay rate F of a metastable state is given in semiclassical theory by 



r = Aexp(-5//l)[l + 0{h)] , (30) 

where A and B depend on the specific shape of the potential. The investigation carried out 
so far allows us to identify the coefficients A and B in Eq. (I5U]) with h\J | Det [/i] / Det[0] | / L 
and A[xc/] respectively. Then, the general expression for the finite time/temperature T{T*) 
is: 



[Xcl] 

(31) 

Eq. ( 13T1) is plotted in Fig. [5] against temperature up to T* for three oscillator energies. 
While at low T*, T(T*) merges with the constant decay rate of the infinite time theory, 
an increase of r(T*) is found in all plots above T* ~ T*/2 where the combined effects 
of quantum fluctuations and classical action softening become evident. Approaching the 
crossover, T{T*) deviates from the T = result and reaches a peak value T{Tp) which is 
larger for lower uj. This effect is mainly ascribable to the soft eigenvalue Si. Note however 
that for huj = lOmeV, T(Tp)/f\w ~ 1, signalling that the application of the semiclassical 
method itself becomes questionable. In fact, as noted below Eq. (|T5l) . the latter works when 
^[^^d] > ^ and such condition starts to be well fulfllled by the case huj = 20meV as shown 
in Fig. [3l Clearly the u values making the semiclassical method feasible also depend on M 
which has been assumed light in the present discussion. Heavvier particle masses favor the 
condition T{Tp)/hw < 1 and sustain the applicability of the semiclassical method over a 
broader range of u. 

Above Tp the decay rate smoothly merges with the classical Arrhenius factor as 
A[xci]/h — * V{a)/ KbT* . The temperatures T|, corresponding to the symbols in Fig. [5l 
mark the effective values at which quantum and thermal decay rates overlap. Beyond 
and approaching T*, the decay rate falls to zero as T{T*) oc {T* — T*)^/"^ and the quantum 
tunneling ceases to exist. The latter power law dependence is driven by \J N^"^ / Det^[C)\ oc p. 
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r(r* 



h 



MN- 

27Th 



A 



Det[h] 
Det^p] 



The increase found for the quantum decay rate up to Tp and the subsequent sharp drop 
is interesting also in view of a comparison with activated systems described by classical 



Ginzburg-Landau finite size models 
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33| in which spatio-temporal noise induces transi- 



tions between locally stable states of a nonlinear potential 



34| . The changes in radius and 



;he stability conditions of metastable metallic nanowires are an example of current interest 



35l . l36l | . In classical systems of finite size L a power-law divergence in the escape rate (with 
critical exponent 1/2) is predicted once a critical lenghtscale Lc is approached at fixed T js^]- 
Instead, the quantum decay rate of Eq. (I3T1) cannot be divergent as the small parameter is 
h which, unlike the noise in classical systems, cannot be varied as a function of L at fixed T 



or vice- versa 



]. Accordingly the quantum tunneling decay rate is small and continuous. 
Finally, it is worth pointing out that the decay rate may be computed independently of the 
squared norm A^~^ as the latter cancels out in Eq. fl3Tl) by explicitly inserting Eq. fl22|) . For 
this reason the behavior of the decay rate essentially depends on Det^[0]/ N^'^ consistently 
with the quadratic approximation for the quantum fluctuations which enters the calculation 
at two stages: a) ii determines the form of the quantum action in Eq. (JT6l) and accordingly 
leads to Eq. ( I20i) : h) it allows us to replace the inverse zero mode eigenvalue by the squared 
norm of the bounce velocity, via Eq. fl2T|) . 



7. Conclusion 



I have developed the finite time (temperature) semiclassical theory for the quantum 
decay rate of a particle in the metastable state of a cubic potential model. In the Euclidean 
path integral formalism, the optimal escape trajectory emerges as the solution of the Euler- 
Lagrange equation in terms of Jacobian elliptic functions. Such solution is a time dependent 
bounce whose periodicity naturally leads to relate the temperature T* to the energy of the 
classical motion. Consistently one defines the crossover temperature T* between quantum 
and activated regimes which depends only on the fundamental oscillator frequency uj. As the 
path integral has been solved treating the quantum fluctuations in quadratic approximation, 
the calculations are confined to the low uj range, that is to the low temperature regime. In 
the numerical analysis I have considered a light particle mass and established, for this case, 
the lowest bound of uo values which make the semiclassical method reliable. The stumbling 
block in the calculation of the quantum decay rate is the estimate of the quantum fluctuation 
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effect in tfie finite time tlieory. In particular, I fiave i) derived a compact expression for 
the overall fluctuation contribution to the path integral in terms of the complete elliptic 
integrals and ii) solved the periodic stability equation which yields the low lying fluctuation 
eigenmodes and eigenvalues in polynomial form. The latter point permits to quantify the 
softening of the lowest positive and of the ground state (in absolute value) eigenvalues as 
T* is approached. The softening of the lowest positive eigenvalue is mainly responsible for 
the enhancement in the quantum decay rate above the prediction of the infinite time (zero 
temperature) theory. The behavior of the decay rate has been studied in detail below T*. 
At T* ~ T*, the thermal activation sets in while the quantum decay rate drops to zero 
according to the power law T(T*) oc [T* — T*)^/"^. Similar conclusions may be drawn by the 
analyses of a quartic metastable potential although the decay rate of the latter is smaller 
than in a cubic potential having the same structural parameters. 
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FIG. 1: (Color online) (a) Cubic potential in the real time representation, (b) Cubic potential in 
the imaginary time representation, uj is in meV. The intersections with the constant energy E (in 
meV) define the turning points for the classical motion. 
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FIG. 3: (Color online) (a) Classical Action (in units h) versus Temperature; (b) Squared Norm of 
the bounce velocity times particle Mass versus Temperature, uj is in meV. 
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FIG. 5: (Color online) Decay rates (in meV) for a particle in metastable cubic potentials (for three 
oscillator energies in meV.) versus temperature up to T*. The symbols mark the temperatures 
in which the quantum decay rates merge with the classical Arrhenius factors. 
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